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The A 1 ]! — X 1 S + band system of carbon monoxide, which has been detected in six highly red- 
shifted galaxies (z = 1.6 — 2.7), is identified as a novel probe method to search for possible vari- 
ations of the proton-electron mass ratio (fj,) on cosmological time scales. Laboratory wavelengths 
of the spectral lines of the A-X (v,0) bands for v = — 9 have been determined at an accuracy of 
AA/A = 1.5 x 10 -7 through VUV Fourier-transform absorption spectroscopy, providing a compre- 
hensive and accurate zero-redshift data set. For the (0,0) and (1,0) bands, two-photon Doppler-free 
laser spectroscopy has been applied at the 3 x 10~ 8 accuracy level, verifying the absorption data. 
Sensitivity coefficients for a varying /j, have been calculated for the CO A-X bands, so that an 
operational method results to search for /i-variation. 

PACS numbers: 33.20.Lg, 14.20.Dh, 06.20.Jr, 98.80.Es 



I. INTRODUCTION 

The search for a variation of the dimensionless proton- 
to-electron mass ratio fj, = m p /m e on cosmological time 
scales can be performed by comparing molecular absorp- 
tion lines in highly redshifted galaxies with the same lines 
measured in the laboratory. Through detection of H2 and 
HD lines in absorbing galaxies upper limits on \i- variation 
have been deduced, resulting in a A/i//i < 1 x 1CP 5 con- 
straint at redshifts z = 2 — 3.5, corresponding to look- 
back times of 10-12 billion years towards the origin of 
the Universe P, [3] ■ For a reliable comparison a database 
of accurately calibrated laboratory wavelengths must be 
available and this has been accomplished through laser 
spectroscopic investigations of the Lyman and Werner 
band absorption systems in H2 [3, i] and HD @. Fur- 
ther, to make such a comparison operational, sensitivity 
coefficients must be calculated for all lines in the spec- 
trum; these represent the wavelength shift induced on 
a line by a varying fi. Such calculations have been carried 
out for the H2 molecule using different methods [6j-|8| as 
well as for the HD isotopomer [3]. 

Through radio astronomical observations of the NH3 
molecule even tighter constraints at the A/i/ /j, < 1 x 10 -6 
level could be deduced, because the transitions in the 
NH3 molecule exhibit much larger coefficients. How- 
ever, this ammonia method has only been applied in the 
two systems where NH 3 is detected (B02 18+365 [13, EH 
and PKS1830-211 [13]) at redshifts z < 1, or look-back 
times of 6-7 billion years. In addition the methanol 
molecule was identified as a molecule with radio fre- 
quency transitions exhibiting very large sensitivity co- 
efficients due to the internally hindered rotation mode 
in the molecule [13-15]. Such methanol lines were ob- 



served in the single object PKS1830-211 fT^ Jl7| and a 
constraint of Afi/fi < 1 x 10~ 7 was derived [lJ|. Thusfar 
the radio astronomical observations have been limited to 



z<l. 

Because the number of useful H2 high redshift absorber 
systems is less than ten, additional methods are explored 
for constraining fi variation at redshifts z > 1. Recently 
a number of high redshift observations were reported on 
the A x n - X 1 !]" 1 " vacuum ultraviolet absorption system of 
carbon monoxide, first in Q1439+113 at z a b s = 2.42 [l9| . 
then in Q1604+220 at z = 1.64 jH, in Q1237+064 at 

2 = 2.69 21| . and finally in three additional systems 
Q0857+18 at z = 1.73, bl047+205 at z = 1.77 and 
Q1705+354 at z = 2.04 [If]. While these observations, 
all performed with the ESO Very Large Telescope, were 
mainly used to measure the local cosmic background tem- 
perature we propose to use the high resolution spectral 
observations of CO A-X to search for Afj,/ fj, at these red- 
shifts. We note that, in addition, spectra of the CO A- 
X system as observed toward Gamma-ray bursts could 
be used for the same purpose, although the only system 
with these spectral features detected so far (GRB060807) 
was observed at a too limited spectral resolution [23j j . 
The rest-frame wavelengths of the CO A-X bands are 
in the wavelength range 130 — 154 nm, hence longward 
of Lyman-a, so that the CO spectral features in typ- 
ical quasar spectra will fall outside the region of the 
Lyman-a-forest (provided that the emission redshift of 
the quasar z em is not too far from the redshift z a bs of the 
intervening galaxy exhibiting the molecular absorption). 
The occurrence of the Lyman-forest lines is a major ob- 
stacle in the search for \i- variation via H2 lines [J 0] ■ 

The present study provides the ingredients to make 
the CO A-X system operational for detecting fi varia- 
tion from quasar absorption spectra. In order to extract 
bounds on Afi/ fj, at the competitive level of < 10~ 5 a lab- 
oratory wavelength data set at an accuracy of AA/A = 

3 x 10~ 7 is required. The A : n - X 1 E + absorption system 
has been investigated over decades, with an important 
comprehensive study by Field et ai. [24| , resulting in the 
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FIG. 1: Recording of the jet absorption spectrum of C O 
with the Fourier-transform vacuum ultraviolet spectrometer 
at the SOLEIL synchrotron, (a) Overview spectra; (b) Detail 
spectrum of the A-X (0,0) band. The line indicated by an (*) 
is a Xe absorption line used for absolute calibration. 

most accurate data by Le Floch and coworkers [25|; the 
data were compiled for astronomical use by Morton and 
Noreau [26| . The accuracy of those data was however lim- 
ited to 0.06 cm" 1 or AA/A = 1 x 10" 6 . This is insufficient 
for an accurate constraint on A/i//i, since the laboratory 
values would make up an essential part of the error bud- 
get. Moreover Drabbels et al. [27J had performed an in- 
tricate multi-step excitation study to derive level energies 
of four (e)-parity levels of A 1 !!, v = 0, J = 1 — 4 at an ac- 
curacy of 0.002 cm" 1 ; these results were found to deviate 
by -0.04 cm" 1 from the classical data [25|, [26j]. This has 
led us to perform two independent studies, a one-photon 
absorption study employing synchrotron radiation and a 
two-photon Doppler-free laser-based excitation study, to 
match the accuracy requirement on the wavelengths for 
the relevant (v — — 9, 0) bands of the CO A-X system. 
In addition a calculation of sensitivity coefficients for 
the spectral lines in the CO A-X system was performed. 



II. FOURIER TRANSFORM SPECTROSCOPY 

The unique vacuum-ultraviolet (VUV) Fourier- 
transform spectrometer (FTS) of the VUV DESIRS 
beamline [28} at the SOLEIL synchrotron, which is 
based on wavefront division interferometry 29], was 
employed to record high-resolution absorption spectra 
of CO in the range 130 — 160 nm. The setup was op- 
timized for obtaining spectra at the highest resolution: 
a molecular slit jet was used to reduce the Doppler 
width, and the FTS was set to its lowest linewidth 
of 0.075 cm -1 , corresponding to a resolving power of 
about 900 000. In Fig. Q] an overview spectrum, as 



well as a detail spectrum is presented showing sharp 
lines in the A 1 n-X 1 E+ (0,0) band at a linewidth of 
0.09 cm" 1 . Wavelength calibration of the FTS spectrum 
is derived from the scan-controlling He-Ne laser [29[ 
and further optimized for the absolute scale referencing 
to a Xenon line at 68 045.156 (3) cm -1 13011 . which is 
based on accurate relative measurements [31( and the 
measurement of a VUV anchor line [Hj]. Based on 
these calibration procedures the uncertainty in the line 
positions of most CO resonances is estimated to be 
within 0.01 cm -1 , corresponding to AA/A = 1.5 x 10" 7 . 
Resulting transition frequencies for the CO A-X (v, 0) 
bands arc listed in Tabic [I] When comparing the results 
from this high-resolution VUV-FTS study with the 
classical spectral data [25|, [26| we find an offset of -0.03 
cm" 1 between the data sets. These offsets are found to 
be in agreement with those of Drabbels et al. [27| for 
the four levels studied. 



III. UV LASER SPECTROSCOPY 

In addition two-photon excitation studies on the CO 
A 1 !! -X 1 S + system were performed employing a pulsed 
dye amplifier (PDA)-laser, injection-seeded by the out- 
put of a continuous- wave ring dye laser, delivering nar- 
rowband output at ~ 300 nm |33l] . Simultaneous record- 
ing of I2 hyperfine lines [34l |35| and transmission peaks 
from a stabilized etalon using the fundamental wave- 
length provides frequency calibration. The excitation 
was performed in a counter-propagating beam geome- 
try to suppress Doppler effects [36|]. In a 2+1' resonance 
enhance multi-photon ionization (REMPI) scheme a sec- 
ond laser at 207 nm is employed for ionizing the A 1 II 
excited state population; this laser is pulse delayed by 
~ 10 ns with respect to the first laser pulse to avoid AC- 
Stark effects by the ionizing laser pulse. Delayed pulsed- 
electric-field ion extraction is also employed to minimize 
DC-Stark effects. Figure shows a recording of an ex- 
citation spectrum of part of the A-X (0,0) band. Due 
to availability and wavelength tunability of the ring dye 
laser the laser studies were performed on the A-X (1,0) 
and (0,0) bands. 

The error budget in Table [II] lists the estimate of un- 
certainty contributions from various sources. AC-Stark 
shift corrections were obtained from intensity-dependent 
measurements and extrapolating the transition frequen- 
cies to zero intensity levels. The accuracy of the obtained 
transition frequencies from the laser measurements are 
limited by the frequency chirp in the PDA-laser. Based 
on previous characterization of the chirp of our PDA sys- 
tem [33[, we estimate a 60 MHz contribution from the 
chirp for this study since a lower-order frequency upcon- 
version of the laser was used. Furthermore, the mea- 
surements were performed at lower PDA-pump energies 
compared to Ref. [Hj], which also reduce the chirp effect. 
We estimate the uncertainty of the line frequencies from 
the PDA investigation to be 70 MHz or 0.002 cm" 1 . 
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TABLE I: Calibrated 1-photon transition frequencies (in vacuum cm -1 ) in the A IT -X E + (t>,0) bands. Lines indicated with 
(*) are derived from the Doppler-free 2-photon laser excitation study of the (0,0) and (1,0) bands and ground state level 
energies [33]. The other lines are derived from the VUV-FTS study. 
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data show agreement between both data sets at a statis- 
tical uncertainty within 0.01 cm . Hence, we conclude 
that the calibration procedures on the VUV-FTS spectra 
yield transition frequencies within the estimated uncer- 
tainty of 0.01 cm -1 . This provides a database of CO 
A-X lines for bands (v = 2 - 8) at AA/A = 1.5 x 10~ 7 
accuracy. A slightly worse uncertainty estimate of 0.02 
cm -1 is quoted for the (9,0) band which has worse SNR 
compared to the other bands due to its lower transition 
strengths. 



IV. SENSITIVITY COEFFICIENTS 

For extracting a possible variation of the proton-to- 
electron mass ratio fi a set of highly redshifted wave- 



In view of parity selection rules in the two-photon laser 
experiment opposite parity A-doublct components in the 
A 1 !! state are excited, when compared to the one-photon 
absorption experiment. Accurate values for the A 1 II A- 
doublet splittings can be deduced from measurements 
in the S, Q, R and P branches (the O branch was not 
recorded), and from the accurately known ground state 
level energies [37|. Based on this analysis the laser data 
were converted to transition frequencies for one-photon 
A-X bands at an accuracy of 0.002 cm -1 , and listed in 
Table H 

The results from the laser-based study serve a two-fold 
purpose. First, they provide a reliable database for com- 
parison with quasar data at the level of AA/A = 3 x 10~ 8 
for the A-X (0,0) and (1,0) bands. Second, a compar- 
ison between the present laser data and the VUV-FTS 
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FIG. 2: Recording of part of the two-photon laser excitation 
spectrum of the A 1 !! - X 1 ^ (0,0) band of CO. The etalon 
markers and the I2 hyperfine line (*) recorded in the funda- 
mental are used for frequency calibration. 



TABLE II: Estimated uncertainties (in MHz) for the transi- 
tion frequencies measured with the PDA-laser system. 



Source 



Uncertainty (MHz) 



Line-fitting 
Residual Doppler 
I2 calibration 
etalon nonlinearity 
AC-Stark shift 
DC-Stark shift 
PDA chirp 



2 

< 1 
2 
2 

20 

< 1 
60 



Statistical 


30 


Total 


70 



lengths A* are compared with a set of rest-frame wave- 
lengths Aq via 0: 



% = {l + z abs ){l + K^) 



(1) 



where z a {, s is the redshift of the intervening galaxy ab- 
sorbing CO and the sensitivity coefficient for each 
line to a variation of \x: 



K % = 



din X 1 
iiln/i 



El 



fj { dEj dE* 
— E l g d/j dfi 



(2) 



with El and E l g the excited and ground state energies 
connecting an optical transition in the molecule. It is 
noted that the definition of Eq. ([2]) yields coefficients 
of opposite sign with respect to other studies, where the 
definition Avjv = K^A^l/^l is used (see e.g. Ref. [U, 
UH). The specific definition of Eq. @ is consistent with 
its use in Eq. ([IJ . The relative variation of the proton- 
electron mass ratio is defined as A/z = /i 2 — (1q , meaning 
that a positive value for A/i/ (j, indicates a larger value of 
/j, in the cosmological past. 



Rovibrational level energies E(v , J) can be expressed 
in terms of a Dunham expansion 



E(v, J) = £ Y k ,i (v+l) [J(J + I) A 2 ] 



h-.i 



(3) 



with Yk,i the Dunham coefficients, known to sufficient ac- 
curacy for the X 1 E + state [38| and for the A 1 II state [39j . 
A represents the electronic angular momentum, A = for 
X 1 !^ and A = I for A 1 II. The advantage of the Dunham 
representation of molecular states is that the coefficients 
scale as Yk t i oc [i r \ d fc ^ 2 , with /i re d the reduced mass of the 
molecule. In studies focusing on /j,- variation it is assumed 
that all atomic masses scale as the proton, i.e protons 
and neutrons treated equally, hence jj, re d can be replaced 
by fi in this analysis [1, |4(| ■ Hence the derivatives for the 
level energies can be analytically taken with 



dYkj 
df 1 



Yk, 



, k 
I + - 

[i \ 2 



(4) 



Substitution in Eqs. ([3]) and ([2]) then straightforwardly 
yields the sensitivity coefficients for the A-X lines of 
CO. 

Interactions of the A 1 II state with triplet states in the 
CO molecule perturb the level structure 24]. The level 
shifts play an important role in the comparative analy- 
sis of quasar data to extract fi variation, but these are 
implicitly included in the experimental determination of 
spectral line positions, i.e. in the values of Table HJ How- 
ever, in the calculation of coefficients the admixtures 
of perturbing character into the wave function compo- 
sition in A 1 n states should also be accounted for. In 
Refs. [1, |4(| a model is proposed to calculate the 
coefficients in the case of non-adiabatic mixing between 
B 1 E+ and C 1 n u levels in H2. This model can be adopted 
for CO and in good approximation one may derive coef- 
ficients 



E 



(Xpert Kpert 



(•5) 



where K pure refers to coefficients for transitions of 
the state in consideration and K pert refer to those of the 
perturber states, while a refers to the admixture in the 
wave function composition (given in Ref. [26|). 

Perturbation-corrected coefficients were calculated 
using Dunham representations for the perturbing d 3 A 
state from Ref. [4l|, and e 3 E~, I 1 ^ - , a /3 E + states from 
Ref. [42| . (we have collected the Dunham coefficients 
in the Supplementary material 43]). The second-order 
mass-dependent effect (adiabatic correction) is estimated 
to be at the 5 x 10~ 5 level and has been neglected. 

As follows from Eq. ([5]) the correction to the coef- 
ficients depends on both the amount of admixture of the 
perturbing state in the wave function, as well as on the 
value of the coefficients for the pure perturber state. 
The most pronounced perturbation occurs between A^^n 
(v = 1) and d 3 A (v = 5) with the highest triplet ad- 
mixture of 17% for the J' — 1 state. In addition, the 
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TABLE III: sensitivity coefficients for 12 C 16 A 1 II - X x £+ («'-«") bands. The uncertainty is estimated to be better than 
1%. 
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perturber value is more than four times that of the 
pure A state value, resulting in a correction of 55% to the 
value compared to the case when the perturbation is 
not included. The A 1 !! (v — 0) state is only very slightly 
perturbed by the d 3 A (v = 3,4) states with 0.13% wave 
function admixture. However, the perturber is 10 
times larger than that of the pure A state value, resulting 
in 1.5% correction for K^. The A 1 ]! (v = 4) rotational 
states can have i^ M -corrections up to 5%, while the A 1 ]! 
(v = 6, J = 8) rotational state correction is almost 10%. 
The rest of the bands have if^-corrections that are less 
than 1%. The coefficients are represented graphically 
in Fig. El where the scatter in the (v = 1, 4 and 6) bands 
show the effect of perturbations. The resulting coef- 
ficients are listed in Table Hill The accuracies in these 
coefficients are estimated to be better than 1%, with the 
dominant uncertainty contributions from the perturba- 
tions (non-adiabatic corrections) which depends on the 
quantum numbers v' ', J' of the upper state. 
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FIG. 3: Calculated coefficients for the spectral lines in the 
A-X system of CO. The spread in values, shown enlarged 
in the insets for the (1,0) and (4,0) bands, illustrate the effect 
of perturber states. 
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V. DISCUSSIONS 

coefficients for CO A-X are thus found in the range 
of values -0.002 to +0.066 giving a spread equal to that 
of the Lyman and Werner bands of H 2 . This makes 
the A-X system of CO a good search ground for putting 
constraints on a possible variation of /i, in particular at 
look-back times of 9.5-11.5 billion years where these fea- 
tures have been detected so far [19l - l22j ] . A set of zero- 
rest frame wavelengths of A-X (v = — 9,0) bands has 
been measured with most transitions at an accuracy of 
AA/A = 1.5 x lO" 7 , and for the (0,0) and (1,0) bands 
even as accurate as 3 x 10~ 8 . The combined data sets 
of rest-frame wavelengths and sensitivity coefficients are 
the basis for a /i-variation analysis on CO spectra from 
high-redshift galaxies, where the analysis will be solely 
constrained by the quality of the astronomical data. So 
far six high redshift galaxies have been observed with 
characteristic CO A-X features. The best example is that 
of a CO spectrum observed at z — 2.69 with a column 
density of logV(CO)= 14.17 cm~ 2 in the sight-line of 
the Q1237+064 quasar system, showing a signal-to-noise 
ratio of 10-40 (over the wavelength range) and resolution 
R ~ 50000 after 8 hrs observation at the ESO Very Large 
Telescope [2l[ . Extended observations with attached Th- 
Ar calibration of the astronomical exposures (not avail- 
able yet for the discovery spectrum of Ref. [21j) would 
result in a competitive result on A/1//1 from the CO 
spectra, i.e. an estimated constraint at < 10~ 5 . 

The objects where CO is detected exhibit large column 
densities of H2 and HD molecules, since extragalactic 
abundance ratios are generally iV(H2)/-/V(CO) > 10 4 and 
V(HD)/7V(CO) ~ 1. Hence, future comprehensive \x 
constraining analyses can be performed from the features 
of all three molecules contained in the same quasar spec- 
trum. In view of the increased number of spectral lines, 



the wider wavelength coverage, and the CO lines lying 
outside the Lyman-a forest, well-calibrated (Th-Ar at- 
tached) and good signal-to-noise observations (SNR 40) 
should lead to constraints at A/i//i < 3 x 10~ 6 . The sys- 
tem Q1237+064 would be the first target of choice. 



VI. CONCLUSIONS 

We have identified the A 1 H—X 1 T, + band system of car- 
bon monoxide a novel probe system to search for possible 
variations of the proton-electron mass ratio (/^) on cosmo- 
logical time scales. Laboratory wavelengths of the spec- 
tral lines of the A-X (v,0) bands for v — — 9 have been 
determined at an accuracy of AA/A = 1.5 x 10~ 7 through 
VUV Fourier-transform absorption spectroscopy, provid- 
ing a comprehensive and accurate zero- redshift data set. 
Two-photon Doppler-free laser spectroscopy has been ap- 
plied for the (0,0) and (1,0) bands, achieving 3 x 10~ 8 
accuracy level. Accurate sensitivity coefficients for a 
varying /i have been calculated for the CO A-X bands, 
where the effect of perturbations have been accounted 
for. It is expected that future fi constraining analyses 
that include H2 and HD and CO transitions from the 
same absorber should result in a more accurate and ro- 
bust constraint for A/1//1. 
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Table I. Dunham coefficients Y k i for the states of 12 C 16 involved in the present analysis: the X 1 S + ground state, the A 1 !! 
excited state and the perturber states. The constants were obtained from [1], [2], [3], and [4] respectively. All values in cm -1 . 
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